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Abstract 

We consider m spinless Fermions in I > m degenerate single- 
particle levels interacting via a /c-body random interaction with Gaus- 
sian probability distribution and k < m in the limit I — > oo (the em- 
bedded A;-body random ensembles). We address the cases of orthog- 
onal and unitary symmetry. We derive a novel eigenvalue expansion 
for the second moment of the Hilbert-space matrix elements of these 
ensembles. Using properties of the expansion and the supersymmetry 
technique, we show that for 2k > m, the average spectrum has the 
shape of a semicircle, and the spectral fluctuations are of Wigner- 
Dyson type. Using a generalization of the binary correlation approx- 
imation, we show that for k <C m <C I, the spectral fluctuations are 
Poissonian. This is consistent with the case k = 1 which can be solved 
explicitly. We construct limiting ensembles which are either fully in- 
tegrable or fully chaotic and show that the fc-body random ensembles 
lie between these two extremes. Combining all these results we find 
that the spectral correlations for the embedded ensembles gradually 
change from Wigner-Dyson for 2k > m to Poissonian for fc<m<I. 

PACS numbers: 02.50.Ey, 05.45. +b, 21.10.-k, 24.60.Lz, 72.80.Ng 
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1 Introduction 



The stochastic behavior displayed by spectra and wave functions of quan- 
tum many-body systems (atoms, molecules, atomic nuclei, quantum dots) is 
usually modelled in terms of canonical random-matrix theory (RMT); the 
application of RMT to such systems has been very successful in many cases 
]T|, Q. At the same time, this type of modelling is not completely realistic: 
All many-body systems that occur in nature are effectively governed by one- 
and two-body forces, while the use of canonical RMT is tantamount to as- 
suming many-body forces between the constituents. Indeed, assuming that 
the one- and two-body interaction can be modelled stochastically, one finds 
that the resulting number of independent random variables is much smaller 
than in RMT. For instance, the number of independent two-body matrix 
elements in a shell-model calculation in atoms or nuclei is typically much 
smaller than the dimension N of the matrices involved, while the number 
of independent random variables in RMT is of order N 2 . This poses the 
question whether a more realistic stochastic modelling of many-body sys- 
tems might yield results which differ from RMT predictions. The question 
was addressed in the 1970's with the help of numerical simulations using 
matrices of fairly small dimensions. The main results were: In a certain 
limit, the average level density does not have the shape of a semicircle but is 
Gaussian; the ensembles are neither stationary nor ergodic; unfolding of the 
spectra yields Wigner-Dyson spectral fluctuation properties, see Refs. to 
fH . Interest in model Hamiltonians with random two-body interactions has 
resurged in recent years in several areas of many-body physics (see Refs. (10 



to 1^3), and the question of possible further differences between such models 
and RMT has resurfaced. 

The paradigmatic models for this type of question are the /c-body em- 
bedded Gaussian random ensembles EGOE(/c) and EGUE(fc) introduced by 
Mon and French in 1975 [^8 and fully defined below. The integer k takes any 
value between 1 and m, the total number of particles. A number of useful 
and interesting results for these ensembles have been obtained p8| , [L], |29|] . In 
particular, it was shown that in the dilute limit (number m of Fermions small 
in comparison with the number I of single-particle states) and for k <C m, 
the shape of the average spectrum is Gaussian (rather than a semicircle as 
in the case of RMT). However, a comprehensive analytical approach to the 
spectral fluctuation properties of these ensembles is still lacking. To quote 
from Ref. Pf: "... there is no real theory yet for EGOE fluctuations, the gap 
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being one of the most significant ones in the entire subject". As far as we 
know, the situation has not changed in the last twenty years. The gap is due 
to the fact that, in contrast to RMT, the embedded ensembles do not possess 
the (unitary, orthogonal or symplectic) invariance in Hilbert space which is 
so essential for the successful analytical treatment of RMT. 

In this paper, we report on a novel approach to the /c-body embedded 
Gaussian ensembles. Although physical interest focusses on the case k = 2, 
the generic features are better understood by addressing the case of arbi- 
trary k with 1 < k < m. As always in RMT, we consider the limit of 
infinite matrix dimension, realized by letting I — > oo. One essential tool in 
our analysis is provided by the orthogonal or unitary invariance of EGOE(/;;) 
or EGUE(A;) with respect to the single-particle states from which Hilbert 
space is constructed. We show that this invariance takes over the role of the 
(orthogonal or unitary) invariance of standard RMT. We first focus attention 
on the second moment of the random fc-body interaction in the Hilbert 
space of many-body states. We introduce a novel eigenvector expansion for 
A^ k \ This allows us to derive statements about the shape of the average 
spectrum as a function of k and m. Further information on the spectral 
behavior of EGOE(/c) and EGUE(fc) is obtained with the help of the super- 
symmetry method f3(J |3l| which applies thanks to the eigenvector expansion 
of A^ k >. We derive and solve the saddle-point equations for EGOE(fc) and 
EGUE(/c). We show that the mean level density and the spectral fluctua- 
tions coincide with those of canonical RMT. Most importantly, we calculate 
the first non-vanishing term in the loop expansion and thereby estimate the 
range of validity of the saddle-point approximation. This range is given by 
2k > m. This result is consistent with an estimate obtained by calculating 
the second and fourth moments of the interaction. In order to gain access 
to the regime 2k < m, we generalize the "binary correlation" method of 



Mon and French to the calculation of level correlation functions. The 
method applies in the limit k <C m <C I. We show that in this limit, the 
spectral fluctuations are Poissonian. This is consistent with explicit results 
for k = 1. To reach a physical understanding of our results, we construct lim- 
iting ensembles of random matrices which are either fully integrable or fully 
chaotic. We show that EGOE(fc) and EGUE(/c) lie between these two ex- 
tremes. In this way, we show that - contrary to common belief - for 2k < m 
and / — > oo, EGOE(fc) and EGUE(/c) do not follow Wigner-Dyson statistics, 
and we give a simple argument why this is so. We comment on the apparent 
agreement between canonical RMT and data from numerical simulation and 
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from atomic and nuclear physics and argue that this agreement is caused by 
a finite-size effect which disappears in the limit I — > oo. 

The embedded ensembles and the limiting ensembles are defined in Sec- 
tion [| The second moment is introduced, and an important identity 
derived, in Section [| The eigenvector expansion derived in Section |] is a 
central analytical result and enables us to calculate the low moments of the 
interaction (Section ||) and to apply the supersymmetry method (Section ||). 
The case /r -C m < Hs studied in Section ^ while Section |] contains the 
conclusions. A preliminary account of some of the results has appeared in 



Refs. |a and P| . 



2 Definitions and Elementary Facts 
2.1 Definitions 

We consider m spinless Fermions in I > m degenerate single-particle states 
with associated creation and annihilation operators and a, where j = 
1, . . . , I. The ratio / = m/l is called the filling factor. For reasons explained 



in Section |3]3|, we confine ourselves to the case of less than half filling, / < 
1/2 or 2m < I. Hilbert space is spanned by orthonormal vectors labelled 
\u), . . .. Each such vector is given by a Slater determinant . . . a] ro |0) 
where |0) denotes the vacuum. Uniqueness is guaranteed by the requirement 
Ji < J2 < • • • < 3m- The dimension iV of this space is 

*=(')■ a) 



The random fc-body interaction with k < m has the form 

V k = Yl Oi ./;:<: ■ ■ ■ 4^ " " " fl n ' ( 2 ) 

l<il<32<— <Jfc<I 
l<i^<i2<...<i/ c <i 

We refer to k as to the rank of the interaction. As in the canonical case, we 
use the labels /3 — 1 and (3 = 2 for the orthogonal and the unitary ensemble, 
respectively. For (3=1, the matrix elements are real and obey 

Vj 1 ,...,j k ;i 1 ,...,i k = V j 1 ,...,j h -,i 1 ,...,i k = V i 1 ,...,i k ;3i,-d k (P = 1 ) • ( 3 ) 

For (3 = 2, the matrix elements are complex and obey 

v ji,-,jk;n,-,h = v h,...,i k ;ji,...,jk = 2 ) • (^) 
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Matrix elements not linked by the symmetry relations (Q) or (Q) are uncorre- 
cted Gaussian distributed random variables with zero mean and a common 
second moment denoted by v$. Thus, 

+5 01 S jifi ...8 jkfk 8 iA ...6 iA ) . (5) 

The overbar denotes the average over the ensemble. 

We "embed" the fc-body interaction Vk(f3) in the m-particle space by 
taking matrix elements (z/|\4|/z). For /3 — 1 (f3 — 2), this defines the k- 
body embedded Gaussian orthogonal (unitary) ensemble of random matrices, 
respectively, in short EGOE(A;) and EGUE(A;). We note that for k = m, we 
have done no more than rewrite the canonical ensembles GOE and GUE in 
what may appear a complicated way. Naturally, interest focusses on the case 
k < m. However, it is always useful to check that for k = m our results 
coincide with canonical random-matrix theory. 

Eqs. and (|5|) are invariant under arbitrary orthogonal ((3 = 1) 

or unitary ((3 = 2) transformations of the single-particle states and so are, 
therefore, the ensembles EGOE(fc) and EGUE(fc). 

The choice of determines the energy scale of our problem. Without 
loss of generality, we put Vq = 1 in the sequel. 

2.2 Regular Graphs 

A graphical representation of the ensemble is obtained by assigning to each 
Hilbert space vector a vertex /i, and to each non-diagonal matrix element 
(^|Vfc|/x) which is not identically equal to zero, a link connecting the vertices 
fi and v . The diagonal matrix elements are represented by loops 

attached to the vertices fi. It is instructive to explore the properties of the 
resulting topological structure. 

The number of vertices is obviously given by N, the dimension of Hilbert 
space. The number M of links emanating from any given vertex is the same 
for all vertices and is given by 

To see this, we observe that for the non-diagonal matrix element (z/|Vfc|/z) 
not to vanish identically, the vectors and \v) may differ in the occupation 
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of as many as s single-particle states where s = 1,2, ... ,k. The number M 
is obtained by summing over all these possibilities. With the exception of 
the loops, the resulting structure is what is called a "regular graph" in the 



mathematical literature |54j]. When we represent the structure as a matrix, 
the number of non-diagonal matrix elements which do not vanish identically 
is equal in each row and each column and given by M. For k < m, we have 
M < N — 1 while M = N — 1 for k = m. It is also of interest to calculate the 
total number P of independent links. This number is given by the number 
of matrix elements above the main diagonal which do not vanish identically. 
Hence, 

P = -MN. (7) 

The minimum number L min of links needed to get from an arbitrary initial 
vertex /i to an arbitrary final vertex v is bounded by 



m 
J 



(8) 



The bracket denotes the largest integer smaller than or equal to m/k. The 
inequality <M) holds because the maximum number of particles that must be 
moved, is m. The maximum number of particles that can be moved per link 
is k. For k = m we get L min = 1. We observe that the graph is very densely 
connected: For / 3> 1, L min is much smaller than N. Obviously, N, M, P 
and L min do not depend on (3. 

The number Kp of independent random variables is different for the uni- 
tary and the orthogonal ensemble and essentially given by the square of the 
number of independent choices of the indices ji, ■ ■ ■ ,jk subject to the con- 
straint ji < j 2 < ■ ■ ■ < jk- More precisely, 



K 



(5 (I 



/3 



2 \k 



+ 5, 



pi 



(9) 



The ratio Kpj P of the number Kp of independent random variables and the 
number P of different links is, for k <C m, much smaller than one and, for 
/ < 1/2, approaches the value f3(N + 6pi)/(N — 1) monotonically from below 
as k approaches m. This shows that for k -C m there are strong correlations 
between matrix elements on different links. The correlations disappear as we 
approach the canonical limit k = m. Fig. [l] shows log 10 K\jP versus k/m for 
various values of the parameters I and m. 
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Figure 1: The ratio K\/P, see Eqs. (0) and (|9|), on a logarithmic scale versus 
k/m. Panel (a): / = m/l = 1/5. Panel (b): / = m/l = 1/2. Panel (c): 
/ = 100. Panel (d): m = 12. 
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2.3 Limiting Ensembles 

Fig. |I] suggests that EGOE(fc) and EGUE(/c) both lie between two limits, 
obtained by assigning the minimum and maximum number, respectively, of 
independent random variables to the graph structure defined by the oper- 
ators appearing in Eq. (0). It is helpful to study the properties of these 
limiting ensembles. The minimum number of random variables is obviously 
one. The corresponding ensemble EGE m j n (/c) is given by the matrix elements 
(v\V™ m \fi) of the interaction 

Vr = v Yl a] l ...a] k a ik ...a il . (10) 

'•• ** i •.> •- h > 

l<il<i 2 <...<i k <l 

The factor v is a Gaussian random variable and real (complex) for (3 = 1 
{(5 = 2, respectively). Without loss of generality we may, however, put 
| v | 2 = 1, removing the distinction between the unitary and the orthogonal 
cases. The ensembles EGOE max (A;) and EGUE max (fc) containing the max- 
imum number of independent random variables are obtained by assigning 
(within the constraints imposed by symmetry) a different random variable 
v Vil to each link in the graph. The interaction has the matrix elements 

v vli Yl {v\a ] h . . . a) k a ik . . . a h \v) . (11) 

l<jl<J2<- ■<Jfc< i 
l<i 1 <i 2 <-..<i k <l 

For (3 = 1 ((3 = 2), the matrix v Vil is real symmetric (complex Hermitean, 
respectively). Elements not connected by symmetry are uncorrelated Gaus- 
sian random variables with mean value zero and variance v vp v v i^i = 6^5^^ + 

The ensemble EGE m ; n (/c) is fully integrable and has spectral fluctuations 
which are not of Wigner-Dyson type. We demonstrate this for the two ex- 
treme cases k = 1 and k = m. For k = 1, the interaction has the form 
(X)j=i a j)(X]i=i a i) which obviously commutes with all permutations of the 
labels j = 1, . . . , I and i = 1, ... ,1 of the single-particle states. Hence, the 
eigenvectors of the interaction are simultaneously eigenvectors of irreducible 
representations of the symmetric group: The system is integrable, and the 
spectral fluctuations are not Wigner-Dyson. The Gaussian distribution of 
the factor v in Eq. ( |T0"D does not affect this statement as it only smears the 
entire spectrum of a given representation of the ensemble over an interval 
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of width unity. We have not succeeded in diagonalizing the interaction ma- 
trix analytically for v = 1. However, extensive numerical calculations have 
consistently yielded two different degenerate eigenvalues Ai, A2 which are re- 
produced by the formulas Ai = I and A2 = 0, with degeneracies n\ = (iT 1 ,) 
and n 2 = 

For k — m, it is easy to see that the matrix representation of EGE m i n (m) 
in Hilbert space carries the entry v on every element. Diagonalization of this 
matrix is trivial and gives the eigenvalues Nv (non-degenerate) and zero 
((N — l)-fold degenerate). Again, the ensemble is fully integrable, and the 
spectral fluctuations are not Wigner-Dyson. 

We show in Section ^| that the spectral fluctuation properties of the en- 
sembles EGOE max (A;) and EGUE max (fc) coincide with those of the GOE and 
GUE, respectively. (In actual fact, we do so only for EGXJE max (k)) . Given 
this fact, we see that the two limiting ensembles cover the extreme cases of 
fully integrable and fully chaotic systems, respectively. Fig. |l] then suggests 
that as k increases, the spectral fluctuations of both EGOE(fc) and EGUE(A;) 
may undergo a gradual transition from Poissonian to Wigner-Dyson behav- 
ior. In what follows, we present conclusive evidence for this statement. 



3 The Second Moment of the Interaction 

After these preliminaries, we turn to the central object of the theory: The 
second moment of the matrix (u\Vk(P)\fj) . For (3 = 2 we derive an identity 
for this quantity which is very helpful for what follows. 



3.1 The Second Moment 

By virtue of the randomness of Vk((3), the matrix (v\Vk((3)\n) is a random 
matrix with a Gaussian probability distribution and zero mean value. The 
distribution is completely specified in terms of the second moment A$ )f>a -(/3) 
defined by 

= WWMWMW) ■ (12) 

Using Eqs. (@),(H) and v 2 = 1, we find 

4t<^) = J2 [(v\ a i--- a l a ^--- a n\ (T ) 

l<Jl<J2<---<3fc< ! 
l<i 1 <i 2 <---<ifc<' 
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x[(pl4 ■■■4 h a 3 k --- a h\ v ) 



(13) 



The matrix elements appearing on the right-hand side of Eq. (|T^) are man- 
ifestly real. Using this fact and the identity (p\H\u)* = (v\H*\p) valid for 
any operator H, we can rewrite Eq. fll3"l) in the form 



4%M 



(p\a f n . . . a] k a ik . . . a h \a) 



l<Jl<J2<--<Jfc<i 
l<i 1 <i 2 <---<!fe<i 



x[(p|< ■■■ a l k a j k 



a h \u) 



+<WH4 ■■■4 h a 3k ■■■ a h \p)\ 



(14) 



Eq. ([14]) is the central relation for the theory developed in this paper. As 
a check, we consider Eq. (|14]) for k = m and (3 = 2. There is only one 
term in the sum over the indices i a , a = 1, . . . , m for which aj m . . . Oj 1 \a) 
gives a non-zero result. This is the term by which all the m Fermions in 
\a) are annihilated. The same combination of indices also appears in the 
expression (p\a\ i ...a\ . For this expression not to vanish, we must have 
a = p. Likewise, we conclude that = unless p = u, and similarly for 
the last term {(3 = 1) in Eq. (0). As a result, we find 



vu 



(15) 



This is the result of canonical random-matrix theory, usually written in the 
form (\ 2 /N)(5p U 5p a + S^ppSua) where 2A gives the radius of the semicircle. 
We note that with our normalization we have A = yN. 



3.2 Duality 

The duality relation is very helpful later on in obtaining explicit expressions 
for several quantities involving traces of powers of Vk((3), and in the appli- 
cation of the supersymmetry formalism. We consider the case (3 = 2. For 
arbitrary k < m, the duality relation has the form 

4tU( 2 )=4^ ) (2). (16) 

This relation connects the matrix elements of the fc-body interaction with 
those of the (m — /c)-body interaction. The crucial difference is in the se- 
quence of state vectors labelled p, u, p, a. The duality relation is obtained 
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from a generalization of the argument leading, for k = m, from Eq. (13) to 
Eq. (p!5[). We consider a fixed term in the sum over the j's and i's in Eq. ([14]). 
For A < jS tPa (2) not to vanish, the single-particle states labelled ji, ■ ■ ■ ,jk must 
be occupied in both \p) and \v), and the states labelled i\, . . . ,ik must be 
occupied in both \p) and \a). Moreover, \p) and \a) must have (m — k) 
occupied single-particle states in common, and the same holds for |p) and 
\v). Hence, for each pair of sets of labels j\, ■ ■ ■ ,jk and ii, ■ ■ ■ ,ik for which 
A^ k \2) does not vanish, there exists a pair of sets of (m — k) uniquely defined 
labels pi < . . . < Pm-k and q\ < . . . < q m -k such that the four states labelled 
p, a, p, v have the form 



lA*> = 


(-)H • 


■■ a l a l- 


. .at 


-JO) 


W) = 




t t 


. . al 


-JO) 


\p) = 


(-)M • 


t t 


■ a l 


JO) 




i-) f " a i ■ 


■ a l a l- 


■ - a l 


-JO) 



(17) 

The sets of indices {p a } and {q a } must have zero overlap with the sets {i a } 
and {jo}- The phase factors in Eq. fllTD are not known and will not be 
needed. Eq. (|17D implies 

<Ml< • • • «L fe at m _ fe • • • a tl \v) = , 

(P\4i ■ ■ ■ a \ m -u a rm-k ■ ■ ■ a n W) = (-) /p+/ct , 

(/i|at 1 ...at fc a ifc ...o il |<T> = (-)^ , 
(p\a\ x ...a\ k a jk ...a h \v) = (-y»+f» . (18) 

From Eq. (|18|) we see that the product of the matrix elements in the first 
two equations equals that of the matrix elements in the last two equations. 
We sum the resulting equality over the sets {i a } and {j a } and the sets {p a } 
and {q a }, respectively, using the one-to-one correspondence of both pairs of 
sets. This yields 

4tU( 2 ) = Yl (p\a ] n ...a) k a lk ...a n \a) 



1<J1<J2<- 
l<i 1 <i 2 <--<*fc<i 



x <p|<4 ■ ■ ■ a i a j k ■ ■ ■ a nW) 
E 



l<ri<r2<.-.<r r] 
l<ti <t 2 <...<*_ 
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X (P\4i ■ ■ ■ 4n-* a rm-» • • • a n W) ■ ( 19 ) 

This is Eq. 



3.3 Particle— Hole Symmetry 

In shell-model calculations in atomic and nuclear physics, particle-hole sym- 
metry is often used to simplify the calculations and to relate results for more 
than half-filled shells to those with less than half filling. This procedure 
applies as well in the present case but, as we shall see, is actually of little use 
because we confine ourselves to random interactions with fixed rank k. 

Let 1 6) = U l j=i a )\°) denote the "full shell". There exists a one-to-one 
correspondence between the set of state vectors \fj) containing m particles 



defined in Section pTT], and the set |/2) obtained by creating (m — k) holes in 
the full shell, b\ x . . . b] m _ k |0). The fe^'s with fej = aj are the creation operators 
for holes. We may consider the matrix elements of Vk(j3) in the basis of hole 
states \ ft) rather than in that of the particle states \fi). To make this particle- 
hole transformation effective, we must rewrite the interaction Vk(/3) in Eq. (0) 
in terms of the fe's and fe^'s and bring all hole creation operators up front. 
The last step entails a number of commutators of 6's and tf's. As a result, we 
find that an interaction which has rank k in particle representation, will be 
a sum of interactions of ranks 0, 1, . . . , k in hole representation. This is why 
particle-hole symmetry is not immediately useful in the present context. On 
the other hand, the very use of a random interaction of fixed rank k is an 
idealization, and it is, therefore, legitimate to consider such an interaction 
both for particles (2m < I) and for holes (2(1 — m) < /), the two types 
of interaction not being connected by any symmetry relation. That is the 
point of view taken in this paper. It is then obviously sufficient to study the 
particle case, i.e., to confine m to values 2m < I so that / < 1/2. 



4 Eigenvector Expansion of the Second 
Moment 

Progress in understanding the properties of EGOE(/c) and EGUE(/c) rests on 
an eigenvalue expansion of the second moment A^ k \(3). To this end, we note 
that A^ k \f3) defined in Eq. ( |13|) has the symmetry properties A^^ pu (2) = 
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ApJ^ u {2) and A/J iP(7 (2)* = A^ )tr/ ,(2). Considered as a "grand matrix" carry- 
ing the indices M = (p, v) and M' = (p, a), respectively, the matrix A$ )PI7 (2) 
is, therefore, Hermitean. (We use the term grand matrix to avoid confusion 
with the term supermatrix in the supersymmetry approach). Every Her- 
mitean matrix possesses a bilinear expansion in eigenvectors. We construct 
this eigenvector expansion for A^J M ,(2). In somewhat symbolic notation, the 
eigenvectors obey 

^A$„{2)C$) = A^k)c£> . (20) 

M 1 

Here s is a running label and the index a allows for possible degeneracies. 
When written explicitly in terms of the state labels p, u, p, a, the eigenvectors 
Cj£, turn into matrices and Eq. (p(j|) contains a partial trace, 

E A ^A 2 ) c l7 ] = ^ s ( k )c% a) • (2i) 

per 

The eigenvectors can be orthonormalized (again this involves a trace), 

E^Cg? = N5 st 5 ab ■ (22) 

For reasons which will be obvious later, we have chosen the normalization 
constant as iV rather than unity. The eigenvalue expansion of A { ^ pcT (2) has 
the form 

4*U 2 ) = ^ E k(s \ k ) c % a)c % a) ■ (23) 

sa 

We proceed to construct the eigenvalues A^(k) and eigenvectors . 

In searching for the solutions of Eq. (^TJ) , we are guided by two simple ex- 
amples. The first one is the case of a random potential V(x). Let us suppose 
that the potential has a Gaussian distribution with zero mean value and a 
second moment given by V(x)V(y) = T{x, y). We assume that translational 
and rotational invariance hold and that, therefore, T is a function only of 
\x — y\. Then, it is possible and useful to decompose T into invariants with 
respect to the rotational group. This is accomplished with the help of the 
multipole expansion 

H\x-y\) = Y^iLmM)Yl\®MVt\%,%)r • (24) 

L,M 
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We have used polar coordinates, and the Y L)M 's are the standard spherical 
harmonics. The term Y^(Q X , ^ x )(Y^(Q y , $ y ))* is the projector onto 
an invariant subspace. The spherical harmonics are eigenvectors of JF with 
eigenvalue fi and span the invariant subspace. We look for an analogous 
decomposition of A^ M ,, using the fact that EGUE(/c) is invariant under uni- 
tary transformations of the single-particle basis. We expect that M and M' 
play the role of the polar angles (Q x , $ x ) (or (B y , <& y ), respectively), and that 
the projector Y^(Q X , $ a; )(Y^ f (© 2/ , is replaced by a projector onto 
a subspace of Hilbert space which is invariant under unitary transformations 
of the single-particle states 1, . . . , /. 

Our second example is the GUE. The second moment of the GUE inter- 
action has the form of Eq. flT5|) with (3 = 2. The two Kronecker symbols 
S^vSpcr display the unitary invariance of the GUE in Hilbert space. From the 
point of view of our first example, we look at Eq. (|T^) as kin to a multipole 
expansion. The expansion has only a single term, and the product of the two 
Kronecker symbols is the projector onto an invariant subspace. The subspace 
is spanned by eigenfunction(s) of the second moment. There is one eigen- 
function which belongs to the eigenvalue A 2 = N. It is the matrix 5 ap . All 
traceless matrices are likewise eigenfunctions but belong to eigenvalue zero. 

We turn to Eq. fl2T|) . To construct the eigenvectors, we state a fact implied 
by our derivation of the duality relation: For A^ k '(2) not to vanish, the two 
states \a) and \p) may differ in the occupation of at most (m — k) single- 
particle states, and likewise for and \u). Guided by this observation 
and by the example of the GUE, we begin with the simplest case denoted 
by s = 0. It is easy to check that cffl = 5^ is a normalized eigenvector of 
A^{2) with eigenvalue A(°>(Jfe) = (™) ( l ~ m k +k )- For s = 1, it is likewise easy to 
check that for j ^ i, the vectors C^r given by (fi\a^ai\u) obey Eq. (pID with 
eigenvalue A^(k) = ( m ^ 1 ) ( l ~ m ^ k ~ l ) ■ For j = z, an additional construction 
is needed because (p\Y^j=i a ] a j\ l/ ) = m ^fiu so that this particular linear 
combination reduces to a multiple of the eigenvector C^J . We must, in fact, 
construct traceless matrices. For j = 1, . . . , / these are given by (//|aj-Oj — 
(1/0 2~2p a l a p\ u ) ■ They are eigenfunctions of A^^) with eigenvalue A^'(k). 
Only (Z — 1) of them are linearly independent, however. We observe that 
the eigenvalue A^(k) has degeneracy = (Z 2 — 1). The operators ajctj 

defining the eigenvectors Cffl span a linear space which is invariant under 
unitary transformations of the single-particle basis. (Trivially X^=i a ] a j * s 
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also an invariant). 

It remains to find linear combinations of the eigenvectors which 
obey the orthonormality conditions fl22|) . For j ^ i, we consider the matrices 
(/^|(a]aj + ajaj)|z/) and (n\y/— T(a]aj — a\aj)\i>) of two Hermitean operators. 
These matrices are all orthogonal upon each other. Normalization is assured 
by multiplying each matrix with ^/(l/2)/(/ — l)/(m(Z — m)). For j = i, 

the matrices have the form (/i|ajaj — (1/0 S P a p a pl z/ ) = (/^l-^il^)- These 
matrices are automatically orthogonal to the ones with j ^ %. We note that 
for j, % = 1, . . . ,1, the matrix T^- = (l/N)tr(BiBj) has the form (m//)((/ — 
m)/(l — 1))[1 — (1//)M]. Here, 1 is the unit matrix in I dimensions, and the 
/-dimensional matrix M carries the entry (+1) in every element. The matrix 
Tij has the non-degenerate eigenvalue zero (which arises because ^ Bj = 0) 
and the (I — l)-fold degenerate eigenvalue (m/l) (I — m) /(/ — 1) . Transforming 
Bj with the diagonalizing orthogonal matrix and multiplying the result with 
a/[/(/ — l)]/[m(/ — m)] yields (/ — 1) orthonormal eigenvectors. Thus, we 
have explicitly constructed the D^ l > orthonormal eigenvectors to eigenvalue 
A ( - 1 \k). We observe that these do not depend on the rank k of the interaction. 
The rank affects only the value of A^\k). In order to simplify our notation, 
we will continue to denote these orthonormal eigenvectors by C^i where 
a = l,...,L>«. 

It should now be clear how to proceed: The index s labels the rank of 
the operators appearing in the matrices cjfu\ It is easy to check that for 
{ji, . . . ,j s } with ji < . . . < j s different from {ii, . . . , j s } with i\ < . . . < i s 
(so that no two indices are equal), the matrices (// 1 at . . . a] s aj s . . . a^Jz/) are 
eigenvectors of A^(2) with eigenvalue 

" S ) ■ (25) 

These matrices are orthogonal upon all matrices in classes s' ^ s. They can 
be made mutually orthogonal by using the Hermitean combinations of oper- 
ators (a] x . . . Oj s a is . . . a i± +a\ i . . . a| s % s . . . %) and \^l(a i j 1 . . . a] s a is . . . a h — 
a] . . . a\ s aj s . . . ajj. This is the same construction as used for s — 1. When- 
ever at least two indices from the sets {ji < . . . < j s } and {i\ < ... < j s } 
are equal, the same problem as for s = 1 arises. Indeed, given a matrix 
in class s' < s, we can seemingly "lift" it into class s by multiplying the 
defining operator by Ej=i a j a j] ( ' s ~ s '^ even though the operator Y^j=x a \ a i 
can be replaced by m. Therefore, it is necessary to find linear combinations 
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which are orthogonal upon these constructs. It is not necessary to perform 
the construction explicitly. Suffice it to say that the dimension D^ s ' of the 
resulting linear space is given by = 1 and, for s > 1 by 



A 2 / I 



DW =U -{s-l) ' <26) 

This is because the total number of sets {ji < . . . < j s } and {i\ < . . . < j s } 
is ( ) , and because we have to subtract the number of matrices obtained by 
"lifting" from all lower classes. We denote the orthonormal eigenvectors in 
class s by Cjj^ where a = 1, . . . , D^ s \ Again, the eigenvectors do not depend 
upon the rank k of the interaction, only the eigenvalue does. 

The total number of eigenvectors constructed in this fashion is given by 
Y^o^^ = (m) • This is the dimension of the grand matrix A^(2). We 
conclude that we have found a complete orthonormal set of eigenvectors 
obeying Eq. (|22|), and that, therefore, the expansion (^) applies. We recall 
that the rank k of the interaction appears only in the eigenvalues given by 
Eq. fl25|). Inspection of this equation shows that the sum in Eq. (|23|) actually 
terminates at s = (m — k). For k = m we retrieve the GUE where only s = 
occurs, and where A^ ) (m) = ( l ) = N, in keeping with Eq. ( |15D . 

For the orthogonal case, we observe that the matrix A^ k \l) given by 
Eq. (14) is the sum of two terms each of which have the form of A^(2). We 
apply the arguments for A^ k \2) to each of these terms separately and find 
as a result 

= ^ E A(S) ^ P&CM + CifClf] . (27) 

sa 

The eigenvalues and eigenvectors are the same as for the case (3 = 2. We can 
combine Eqs. ( |2"5| ) and (^) into the single equation 

4 k lM = jj E ^ s \ k )\- c t a)c % a) + ^c ( ; P a)c t a) ] ■ (28) 

sa 

Eqs. (|22|) , fl25|) , (|26|) and (0), constitute the central results of this Section. 

Another important identity is obtained by combining the eigenvalue ex- 
pansion for A( k \2), Eq. (0), with the duality relation (0). This yields 

5>M(*)C£0C£0 = EA ( V - k)C$CW . (29) 

sa tb 
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This equation implies further useful relations. For instance, putting p = v 
and fi = a, summing over /i and p, and using the orthonormality relations 
we find 

A (s) (k)D (s) = iVA {0) (m — k) . (30) 



5 Moments of 



We apply the results obtained in Sections [5] and |] and calculate the depen- 
dence on k, m and I of measures of the shape of the average spectrum of 
EGOE(/c) and EGUE(A;). As always in this paper, we do so in the limit 
/ — > oo. Only in this limit do we expect to obtain generic results. This 
limit corresponds to the limit A" — > oo of canonical random-matrix theory. 
As a by-product, we obtain results which allow us to infer how quickly the 
embedded ensembles approach the limit I — > oo. 

We calculate ensemble-averages over products of traces of powers of 
Vk{(3). Traces are always taken in the Hilbert space of Slater determinants 
\p). From Eqs. fl2S| ) and ( |30"D we have 

m—k 



ItrVM) 2 = £ A (s \k)D^ = (i + M A <°)(m - k) (31) 



' N 



N 2 



s=0 



and 



Ur{V k {P)Y = ^\k) + 5 pi O{l/l 2 ) . 



(32) 



Two further averages are given by 

(A(°)(A;)) 2 + 



-^tr(W))*! 



m—k 



(A(°)(fc)) 2 + 



2(1 + 5, 



01 J 



N 2 



^(A^(m- A;)) 2 L> (s) 



+5 pi O{l/l 2 ) 



(33) 



and 



min(k,m—k) 

-trJMPW = 2(A(°)(fc)) 2 + - £ A^(k)M s \m-k)D^ 



s=0 



17 



+5 pi O{l/l 2 ) 



(34) 



Both these results can easily be found for (3 = 2. For j3 = 1, we en- 
counter terms of the form ^2 u C^Cpt a ^ which cannot be worked out with 
the help of the orthonormality relation (|2"2"|). We recall the explicit con- 
struction of the matrices in Section |j. When the indices ji, ■ ■ ■ ,j s and 
i\, . . . , i s of the operator defining the matrix C^t are all different, the Her- 
mitean operators come in pairs of the form (A + A') and i(A — A') where 
A = (at . . . at aj s . . . a^). The matrix C^t containing (A + A^) is real and 
symmetric, the one containing i(A — A') is purely imaginary and antisym- 
metric. Using these symmetry properties and completeness, we find that in 
the sum C^v Cpt < , the two operators A and A^ appear only in the forms 
(A) 2 and (A^) 2 . These terms both vanish because A is nilpotent. We are 
left with the contributions where at least one pair of indices ji, ■ ■ ■ ,j s and 
ii, . . . , i s of the operator defining the matrix coincides. The number of 
such terms in the sum over a is of relative order I /I 2 and, thus, negligible in 
comparison with other terms which we keep. 

We construct three ratios which yield information about the shape of 
the spectrum. The ratio S measures the fluctuations of the center of the 
spectrum in units of the average width of the spectrum. It is given by 



S = ' . (35) 

^r(V k (f3)) 2 

The relative fluctuation of the width of the spectrum is given by 



Ur(V k ((3)) 2 ) 2 ~ Ut<V k (P)) 2 ) 2 
R = ± —L K — 2 '- . (36) 



The kurtosis k is given by 



2 + Q . (37) 



The quantity Q marks the difference in spectral shape between the semicircle 
(« = 2) and the Gaussian (ft = 3). 
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The ratios S and R yield information about finite-size effects in numer- 
ical calculations of EGOE(fc) and EGUE(A;) p6[ . The ratio S is given by 
S(k, m, I) = ((1 + 6fn)/N)(AW(rn - k)/A^(k)). The ratio R has the value 



s=0 



2(1 + ^) « / AW(m-fc) y (s) 



s=0 



The result for i? is consistent with the estimate by French ||. For k > 1, 
we have both 5 — > and -R — > for / — > oo. However, for fixed values of 
and f,Soc l~ 2k oc (lnAQ~ 2fc and R oc / _2fc oc (lnAQ~ 2fc decrease very slowly 
as I increases, much more slowly than in the case of the canonical ensembles 
where S(m,m,l) = 1/N 2 and R(m,m,l) = 2/N 2 . This fact illustrates the 
difficulty in obtaining reliable spectral information numerically. We note 
that -R(0, m, I) — 2 and return to this point in Section [7| in the context of the 
binary correlation approximation. 

From the point of view of the present study, the most interesting infor- 
mation resides in the ratio Q defined in Eq. (|37|) and given by 

We have Q(0,m,/) = 1 and Q(m,m,l) = 1/N 2 , consistent with a Gaus- 
sian and, for iV — > oo, a semicircular shape of the spectrum, respectively. 
Moreover, recalling that / = m/l, we have the following limits: 

lim^ooQ = , for 2k > m with both k and m fixed , 

(m— k\ 

lim^ooQ = , for 2k < m with both and m fixed , 

lim^ooQ = , with both k/m and / fixed , 

lim^ooQ = 1 , with both k and / fixed . (40) 

A graphical representation of lim/^ooQ is given in Figure 0. We note that 
as m increases, the curves shrink toward the ordinate. In the limit m — > oo, 
the value of Q is unity at k/m = and zero everywhere else. 
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0.0 0.5 1.0 

k/ m 

Figure 2: The limit I — > oo of Q(k, m, I) versus k/m for three values of m. 

These results show that in the limit I — > oo and for fixed m, the spectral 
shape is that of a semicircle as long as 2k > m. (In the unphysical limit where 
with I also k tends to infinity, we also obtain a semicircle). The semicircle 
undergoes a smooth transition to Gaussian shape for 2k < m. We have 
identified the critical value of k where the transition of the average spectral 
shape from semicircle to Gaussian begins. Obviously, the kurtosis alone is 
not sufficient proof for this transition: We have not investigated the higher 
moments of the spectrum. The calculation of such moments is very difficult 
at best. In Section ^| we present independent evidence for the conclusion that 
the transition to Gaussian shape begins at 2k = m. 

For k = m, we have Q oc 1/N for I — > oo. For small k, the approach is 
not quite as rapid. For instance, for k = 1, we have Q oc l 3 /[Nm(l — m)]. 



6 Supersymmetry Approach 



Until now, the embedded ensembles have largely resisted attempts to ap- 



ply the supersymmetry method [2S]. This is because these ensembles lack 



the orthogonal (unitary) invariance in Hilbert space characteristic of GOE 
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(GUE, respectively). The eigenvector expansion ( |2lf ) remedies this deficit: 
After averaging over the ensemble, the integrand of the generating functional 
contains an exponential the argument of which is a sum of squares of bilin- 
ear forms in the integration variables. This expression can be linearized via 
the familiar Hubbard-Stratonovich transformation. We use the saddle-point 
approximation, solve for the saddle-point manifold, and consider the lowest- 
order terms in the loop expansion for the one-point and the two-point func- 
tions. The results in Sections |T2] and |0| supersede and correct statements 
made in Section 5 of Ref. [ 3^| . In a last Section, we show that the spectral 
fluctuations of the limiting ensembles EGOE max (A;) and EGUE max (fc) defined 
in Section ^]3| coincide with those of the GOE and GUE, respectively. For 
simplicity, we display only the case f3 = 2. The calculations for (3 = 1 can be 
carried through in full analogy. 



6.1 Saddle Point 



The steps leading to the saddle-point equation are perfectly standard |3j], |3_L 
and are only sketched here. We consider the case of the two-point function 
and obtain the one-point function later by specialization. The integrand of 
the generating functional contains the factor 

exp [ I £ % pa L 1 £^\V k \v)L 1 J^ vpa ] . (41) 

fiupa 

Here, ^ ^ pa stands for the complex integration variables with a = 0, 1 for 
commuting (anticommuting) variables, p — 1,2 for the advanced (retarded) 
Green function, and fi = 1, . . . , N. The diagonal supermatrix L distinguishes 



the retarded and advanced case and is defined in Ref. [51]. Averaging expres- 
sion PH) over the ensemble and using Eq. (|12|)> we obtain a quartic term in the 
\l/'s involving the matrix A^ k \ We use the eigenvector decomposition (^3|). 
This allows us to perform the Hubbard-Stratonovich transformation. For 
each value of the pair (s, a), we introduce a supermatrix of composite 
variables. Integration over the \l/'s yields an integral over the composite vari- 
ables containing the factor exp(— C e g). The effective Lagrangean £ eff has the 
form 



N 



^trg[a^] 2 +tr /t trgln[(E-ieL)^-^AW(A;)a^^-jJ . (42) 
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The A*- S )(&)'s are the positive square roots of the eigenvalues A^ s \k). The 
energy arguments E\ and E 2 of the advanced and the retarded Green function 
are used to define E = (1/2) (E\ +E 2 ) and e = E 2 — E\ while J stands for the 
source terms. The occurrence of the factor iV ^> 1 in the effective Lagrangean 
causes us to use the saddle-point approximation. The wish to make this 
factor explicit was the reason for the normalization chosen in Eq. ([23|). Under 
omission of the small terms J and e oc 1/N, the saddle-point equations read 



a (sa) 



-itr^^S-^A^^ife)^* 6 ^^]-^^)^) • (43) 



N 

tb 



We have used matrix notation in Hilbert space. The saddle-point Eqs. ([43|) 
constitute a set of coupled equations, one for each of the unknown quantities 



<j( sa ). 



To solve the saddle-point equations (|4"3|), we multiply the equation for 
a( sa ') with \( s \k) and sum over all (s,a). We find 



— jytr M 



( J E-X)- 1 ^A (s) (A;)C (sa) jc^ a) (44) 



where X pcr = ^2 sa X^ s \k)a^ 'CpT '• Comparison with Eq. (^) shows that the 
right-hand side of Eq. (f44"l) can be expressed in terms of the grand matrix 
A^ k \ Doing so and using Eq. (0), we obtain 

X P * = E (p\V k \^)[{E - X)-\MV k \°) ■ (45) 

fJLV 

The overline refers exclusively to the two matrix elements of displayed 
explicitly. We use Eq. (^) in the identity G(E) = [E — X]^ 1 for the averaged 
Green function (a one-point function). This yields the generalization of the 
Pastur equation to EGUE(/c). Iteration with respect to X of this equation 
yields a series which can be interpreted in terms of diagrams: Every averaged 
pair of Vfc's is connected by a contraction line. The series differs from the 
series obtained by averaging the Born series for G(E). It lacks all those 
contributions where at least two contraction lines intersect. We return to 
this point in Section |^. Terminating the iteration after a finite number of 
steps (this number is arbitrary), we find that X obeys the GUE saddle- 
point equation X = A^(k)[E - X)" 1 . We first look for a solution X diag 
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which is diagonal in the superindices. For \E\ < 2X^(k) this yields X^ g = 
d^X^ik)^ with 

^-Jm^f-iS^ 2 - (46) 

The ± signs refer to the retarded (G~(E)) and the advanced (G + (E)) case, 
respectively. Inserting this result back into the saddle-point equations (|43"D , 
we find that is identical to while = for all s > 1. For the 
one-point function this shows that, within the range of validity of the saddle- 
point approximation, the semicircle describes the generic form of the average 
spectrum of EGUE(fc). In contrast to the standard GUE case, the radius 
X^(k) of the semicircle tends to infinity for I — > oo. The mean level spacing 
d(k) oc \<®(k)/N tends to zero (m fixed) or remains finite (/ fixed) in the 
same limit. 

In the case of the canonical ensembles, the invariance of the effective La- 
grangean under general pseudounitary transformations implies that the two- 
point function possesses not a single saddle point but a saddle-point mani- 
fold. The same situation prevails here. The manifold is given by T~ 1 <t^T. 
Here is a diagonal supermatrix of dimension four with entries given by 
Eq. ( pp in the usual way. The matrices T parametrize the coset space 
UOSP(l, 1/1, l)/[UOSP(l, 1)®U0SP(1, 1)]. Using this fact and the saddle- 
point values cr (sa ' ) = for s > 1 in the effective Lagrangean, we find that the 
first -order term in e takes the canonical form 

- — tTg(LT^LT) . (47) 

Here d is the average level spacing. The argument also carries through for 
higher-point correlation functions and shows that the spectral fluctuations 
of EGUE(fc) are identical to those of the GUE. 

6.2 Loop Correction to the One— Point Function 

We wish to determine the range of validity of the saddle-point solution in 
the limit I — *■ oo. We do so by using a power-series expansion for the ef- 
fective Lagrangean around the saddle-point and a subsequent expansion of 
the resulting exponential ("loop expansion"). We recall that in the GUE 
case, all terms of the resulting series save the first vanish individually as the 
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dimension N of the GUE matrix tends to infinity. This is true except at 
the edge of the spectrum. In the present case, the calculation of the terms 
in the loop expansion is a formidable task, and we confine ourselves to the 
correction term of lowest non-vanishing order. We investigate the behavior 
of this term for / — > oo. 

For the one-point function, the sign in Eq. is fixed (negative, say), 
the dimension of the supermatrices is two, and we deal with a single energy 
variable E\ only. We write = 5 s qt^ + 5cr^ sa ' and expand the effective 
Lagrangean in powers of 8a^ sa \ For the logarithmic term, we effectively 
generate an expansion in powers of = J^ sa \^ s \k)C^ saS> Sa^ sa \ In this 
expansion, we keep terms up to fourth order. We expand each of the resulting 
expressions up to first order in the source variable J, and the exponential 
in powers of the terms so generated. Details are given in Appendix We 
obtain 

— — _ iVrWtrgJ / (r (Q))4 A (s) {k)A (s) {m _ k) 
[) A<*»(*) \ + N(Y(k)y^ (A(°W 

The squared factor Y(k) is some suitably chosen complex number depending 
on k, m, /, and E. (Y(k))~ 2 compensates for dropping factors Y^ s \k) = 
1 — (t^ ') 2 (k) / A^ ' (k) in the denominators inside the sum in Eq. (f48|). 
These factors Y^i^k) originate from the Gaussian weight factors. Y{k) is 
well defined with |V"(fc)| finite as long as |Y^(/c)| is finite for all s. At the 
center of the spectrum where (r^) 2 = (—1), we have 1 < Y^^k) < 2. (More 
generally, we have < |y^^(fc)| < 2 everywhere in the spectrum.) In this 
case, Y(k) is real and lies within the same bounds. For the loop correction, 
this gives the upper bound Q(k, m, I), and Eq. (filf) simplifies to 

mO)\<^^-(l + Q(k,m,l)) . (49) 

The loop correction is bounded by the ratio Q defined in Eq. (p9[), cf. 
also the relations (fit]). This fact is highly satisfactory: It agrees with 
our evaluation of the kurtosis in Section [| in the sense that non-vanishing 
corrections to the semicircle shape occur for / — > oo whenever 2k < m 
while the semicircle form prevails for 2k > m. In view of the bounds 
1 < |Y( S )(A;)| < 2 on Y^ s \k), the asymptotic behavior implied by the re- 
lations ( |4*0D holds true independently of the estimate leading to the rela- 
tion (|49|) . And this behavior prevails within the major part of the spectrum. 
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With z^(k) = AW(fc) /A(°)(fc) and x = E / (2X ( -°\ k)), this statement fol- 
lows from \Y^(k)\ = + 2zM(l - 2x 2 ) + (z^)) 2 . The modification of the 
shape of the spectrum due to the loop correction is determined by the energy 
dependence of which contains the term y/1 — x 2 . At the edge of the spec- 
trum, Y^(k) approaches zero. So does Y(k) and a more detailed analysis 
would be required to see whether the square-root singularity is asymptoti- 
cally removed by the loop correction. We have not pursued this point. 

We conclude that the average spectrum of EGUE(/c) has semicircle shape 
for 2k > m while corrections to this form arise for 2k < m. These correc- 
tions may be substantial. Thus, the conclusion reached in Section [| gains 
independent and strong support by the present analysis. 



6.3 Loop Correction to the Two— Point Function 

For the two-point function, we again expand around the saddle point. The 
logarithmic term in the effective Lagrangean has the form 



tr„trg In 



E - A (( V 0) - -eTLT- 1 - TJT' 1 - V 5£ (sa) 
2 ^ 

sa 



(50) 



We have absorbed the matrices T and T~ l in the definition of the four-by- 
four matrices 5T^ sa \ Except for s = these matrices have 16 independent 
variables. For s = 0, we have 5a^ = T5PT' 1 where 5P commutes with 
L. This is because the remaining degrees of freedom in 5a^ generate the 
saddle-point manifold. 

We expand the logarithm in powers of and keep terms up to order 
four. We expand the resulting terms in powers of 1/2 e TLT~ l + TJT^ 1 
and keep terms up to order two. We expand the exponential in powers of 
these terms. After integrating over the Sa's and dropping terms that vanish 
as 1/N for — > oo, we write the result as an exponential, see Appendix 0. 
This exponential contains two terms. The first contains the expression 



me , , 

—trgiLT^LT) . (51) 
«i 

Here d\ is the modified average level spacing which results when the last term 
in Eq. fl48|) is included in the average Green function. The index 1 stands 
for loop. This modification is no more than a consistency check. Indeed, the 
graded trace in Eq. (|5lD yields Wigner-Dyson spectral fluctuations. These, 
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however, must be expressed in terms of the correct average level spacing 
which, in the framework of the loop expansion, is precisely d\. 
Interest focusses on the second term. It is proportional to 

R(k, m, I) trg [1-eLTLT- 1 + LTJT^ 1 ]) 2 . (52) 
Vai 2 / 



In a different context it was shown by Kravtsov and Mirlin [[35] that such 
a term produces deviations from Wigner-Dyson spectral fluctuations. The 
strength of these fluctuations is proportional to the ratio R(k, m, I) given in 
Eqs. (|38|). We recall that R vanishes for I — > oo but does so particularly 
slowly for small values of k and fixed /. 

Whereas the lowest non-vanishing loop correction has given us clear evi- 
dence for a change in spectral shape at 2k = m, this is not the case for the 
spectral fluctuation properties. The ratio R vanishes asymptotically more 
slowly than 1/N but does not indicate a break at 2k = m. If the transition 
to Gaussian shape is accompanied by a change in level statistics, we would 
expect that terms like the one in ( |52"D but which do not vanish asymptoti- 
cally, would arise in higher-order terms of the loop expansion. We return to 



this point in Section 7.1 



6.4 Universality of EGUE max (/c) 

To show universal random-matrix behavior for EGOE max (/c), we apply the 
supersymmetry formalism. The steps are quite similar to those of Sec- 



tions pTT], B72L and 5.3, and we confine ourselves to indicating points of dif- 



ference. We define the operator 

4 = £ «!,...<. (53) 

l<h<h<-<h<l 

This operator is obviously different from the one in Section |5] where the same 
symbol was used. After averaging and the Hubbard-Stratonovich transfor- 
mation, we obtain the effective Lagrangean 



-tr^trgcr 2 + tr^trg In 



(E - -J- ieL)^] . (54) 



The quantity a is both, a supermatrix and a matrix in the Hilbert space 
with m — k particles. The saddle-point equation for a can be solved as in 
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Section ^TT] and yields the solution {v\A^ a^ A\v) = Tm2 x which is independent 
of v. The quantity Tmlx has a form similar to that of in Eq. (|46|) and is 
given by A = E/2±iy/N(k) - (£/2) 2 , with N(k) = (1/N) tr^AU^U) a 
constant and N(k) ^> 1. This shows that the average spectrum has semicircle 
shape with radius 2^N(k). The saddle-point manifold is again generated 
by the matrices T. Expanding the effective Lagrangean to first order in e 
leads to a term of the form P7|). This shows that the level fluctuations 
are of Wigner-Dyson type. All of this is quite straightforward. The crucial 
question is: Do the leading terms in the loop expansion vanish for I — > oo? 

We consider the one-point function and write a = a^ + 5a. The loga- 
rithmic term contains the matrix elements (ulA^SaAlfi) of 5a. The steps are 
the same as in Section [T2] but can be carried through independently for each 
value of /i, cf. Appendix [A]. The correction is proportional to (r?) /N(k)) 4 . 
Even in the worst case k — 1, this contribution vanishes like l~ A for I — > oo. 
For the two-point function, we find likewise that the first non-vanishing 
term in the loop expansion behaves like ((™)/(-^V(^)) 4 - We conclude that 
EGUE max (/c) has the same spectral properties as GUE. This is consistent 
with results by Mirlin and Fyodorov [[JIJ who showed that sparse random 
matrices possess Wigner-Dyson statistics as long as the number of non-zero 
matrix elements per row and column exceeds the value unity. 



7 The Case k < m < / 

The results of the previous Sections show that for 2k > m, EGOE(/c) and 
EGUE(A;) behave generically: The average spectrum has semicircle shape, 
and the spectral fluctuations coincide with those of GOE and GUE, respec- 
tively. We now address the question: What happens in the regime 2k < ml 
We expect a transition to Gaussian spectral shape but have little informa- 
tion yet on the spectral fluctuation properties. In this situation, it may be 
useful to consider far removed as possible from the transition point 

2k = m, i.e., to consider k m. For technical reasons, it is then advanta- 
geous to take also m <^il. This is the case dealt with in the present Section. 



As a preparatory step, we consider in Section [TT] the case l = fc<Cm<C 
/. We show that the spectral fluctuations are Poissonian. By calculating 
the inverse participation ratio, we show, furthermore, that the eigenstates 
display localization in Fock space. Both features are non-generic. It is 
conceivable, of course, that the case k = 1 is special, and that already k = 
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2 behaves very differently. In Section T71, we show that this is not the 
case. We use the "binary correlation appoximation" developed by French 
and collaborators |2S|, [L], p7[| . With the help of this approximation, it was 
shown 
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that for k <C m <C I, the average spectrum has Gaussian shape. 
We generalize the approximation in such a way that it also yields information 
on two- and higher n-point correlation functions. We show that the level 
fluctuations are Poissonian for k -C m <C I. 



7.1 The Case k = 1 

To simplify notation we again confine ourselves to the unitary case. The 
extension to the orthogonal case is perfectly straightforward. For m = 1 the 
random one-body operator can be diagonalized for each realization of the 
ensemble. With a = 1, ...,/, we denote the eigenvalues and eigenfunctions by 
e a and ^/ a , respectively. The eigenfunctions can be expanded in our single- 
particle basis, ^ a = cJjO) = X^=i U a jOj\0). For I — > oo, the eigenvalues 
e a obey GUE statistics, the average level density has semicircle shape, and 
the elements of the matrix U a j are complex Gaussian distributed random 
variables with mean value zero and a second moment given by 

U ai U* pj = - SapSij . (55) 

For terms of higher order, we use the results of Ullah and Porter given in 
Eq. (7.16) of Ref. Q. 

For m > 1, we can easily find the eigenvalues and eigenfunctions of this 
ensemble. The eigenfunctions are orthonormal Slater determinants given by 
\Xn) = 4*1 • • • c a m |0) where the a's take values from 1 to / with «!<...< a m . 
The corresponding eigenvalues are 

m 

E n = J2^- (56) 

It is obvious that form > 1, the average level density has Gaussian shape: It 
is the convolution of m semicircular level densities. It is equally obvious that 
for / > m > 1, the spectral statistics of EGUE(l) is Poissonian. Indeed, 
from Eq. (|56|) we see that neighboring eigenvalues E n and E n+ \ may be com- 
posed of quite different combinations of single-particle energies Ej, and both 
level repulsion and stiffness of the single-particle spectrum lose importance. 
This is, to some extent, true already for m = 2. 
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The example is instructive because for m = 2, we are right at the tran- 
sition point 2k = m separating the regime of the semicircle from that of the 



Gaussian, see Sections pT2| and |673|. The correction fl52|) derived in Section p2 
contains the factor R which vanishes asymptotically for I — > oo whereas the 
case m = 2 discussed above leads us to expect that non-vanishing finite cor- 
rections to the Wigner-Dyson fluctuation behavior exist in this limit. This 
fact reconfirms our expectation that higher-order terms in the loop expansion 
yield contributions of the form (trgLTLT^ 1 ) 2 which do not vanish asymp- 
totically. 

We turn to the eigenfunctions and ask whether these display localization 
in Fock space. This would be another strong indication for non-generic 
behavior. To this end, we expand the eigenfunctions in terms of the Slater 



determinants \n) introduced in Section [TT 

N 



\Xn) = ^29np\li) (57) 



and calculate the inverse participation ratio 



N 



V = Nj2\9n,\ 4 - (58) 



For states which are completely mixed, we expect \g n ^\ 2 — 1/N thus |g n/ J 4 = 
1/N 2 and V m 1. A value of V significantly larger than unity signals local- 
ization. 

The eigenvector \x n ) has m occupied single-particle states labelled a p 
with p = 1, . . . , m and «!<...< a m . For the corresponding indices are 
j q . The expansion coefficients g n/Ji = det(U ap j q ) are determinants formed of 
the matrix elements of the unitary transformation U a j introduced above. It 
is straightforward to show that |g n /J 4 can be written as a determinant, 

|^| 4 = det ( J2 VopiVpPfuVZj) ■ ( 59 ) 

/3 — a l5 ...,a m 

In averaging |g n/ J 4 over the ensemble, we keep the leading terms in an ex- 
pansion in inverse powers of m. The average is calculated using Eq. ( |55l) and 
Wick contraction. The leading contribution in powers of 1/m is obtained by 
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contracting pairs of Vs in such a way that the minimum number of sum- 
mations over the indices % and (3 is destroyed by the Kronecker delta's on 
the right-hand side of Eq. fl55|). Consider an arbitrary entry in the matrix 
appearing in the determinant in Eq. (0). The indices (3 and i of the sec- 
ond factor UL are both summed over, suggesting that this factor should be 
contracted with either of its neighbors to the right and left. Doing so and 
using Eq. (|55D we find for the entry the value 2{m/l) J2i^a p iU*, The fac- 
tor two accounts for both ways of performing the contraction. Ullah-Porter 
corrections to this result are small of order m/l. Hence, 



W 4 ^(— j det£ U avi U* alpi . (60) 

Inserting this result back into the expression (^) for V and using normaliza- 
tion, we obtain V ~ iV(^p) m . We take the limit / — ► oo of V for fixed m » 1 
and find V — > (2e) m ^> 1. We recall that we expect V ~ 1 for delocalized 
states and conclude that for < rn < I, the eigenstates of EGUE(&;) show 
localization in Fock space. This statement applies uniformly to all eigenstates 
and is not connected with a derealization transition. 



7.2 Binary Correlation Approximation 

The results of Section [TT] raise the question: are these results specific for the 
case k = 1 (which then would play a somewhat singular role), or are they 
typical for I; C m C i? With the help of a suitable generalization of the 
binary correlation approximation [^, |l|, p7| , we now show that the second 
alternative holds. We proof that the eigenvalues of the EGUE(A;) have a 
Poissonian distribution in the limit I; < m ■€ 1 in general. 

The binary correlation method has been successfully applied to calculate 
the average one-point function of the EGUE(/c) p8|, In Ref. |29|, the aver- 



age one-point function was expanded in a power series in Vk- One makes use 
of the fact that Vk is a Gaussian distributed random operator and performs 
Wick contractions to carry out the average. In the limit k <C m <C I, each 
pair of Wick-contracted VVs may be approximately replaced by vlA^(k). 
The resulting expression can be resummed using Borel summation. This 



yields the Gaussian spectral shape first found by Mon and French [£8| . 

We proceed similarly for the two-point function by expanding both Green 
functions in powers of Vk and applying Wick contraction. The disconnected 
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part is calculated in exact analogy to Ref. [25]. However, cross-contracted 
pairs of V^s which arise in the connected part have to receive special treat- 
ment. In Appendix [BJ, we show that in the terms of order (Vk) 2n , the contribu- 
tions of cross-contracted pairs carry lower powers of I than the contributions 
to the disconnected part. In the limit k <^ m <^ I, this fact causes the con- 
nected part to vanish relative to the disconnected part. More precisely, the 
properly normalized two-point correlation function 



R^z 2 ) = M?kL-i (6i) 

g(zi) ■ g{z 2 ) 

vanishes in the limit I — > oo. More generally, all higher n-point functions 
turn out to be asymptotically given by the product of n average one-point 
functions. Hence, all connected n-point correlation functions vanish in the 
limit k -C m <C I, yielding Poissonian spectral fluctuations. 

The binary correlation approximation applies for non-zero values of k 
as long as I; < m < (. The approximation predicts Gaussian spectral 
shape and, perhaps unexpectedly, Poissonian spectral fluctuations. This is 
in keeping with the results trivially expected for k = and with the results 
for k = 1 found in Section |7.1] . In terms of spectral properties, the cases 
k = and k = 1 are seen not to play a special role. In the limit k m ^ I, 
the spectral properties do not show discontinuous changes when k approaches 
zero. This corresponds to the continuous dependence on k in all the formulae 
derived in earlier Sections. 



8 Conclusions 

We have studied the shape of the average spectrum and the eigenvalue fluc- 
tuations of the embedded ensembles EGOE(/c) and EGUE(A;) in the limit 
of infinite matrix dimension, attained by letting the number I of degener- 
ate single-particle states go to infinity. We have shown that for sufficiently 
high rank k of the random interaction (2k > m where m is the number of 
Fermions), these ensembles behave generically: The spectrum has semicir- 
cle shape, and the eigenvalue fluctuations obey Wigner-Dyson statistics. A 
smooth transition to a different regime takes place at or near 2k = m. It has 
long been known that the average spectrum changes into Gaussian shape, al- 
though the point of departure from the semicircle shape was not known. We 
have shown that in addition — and contrary to general expectations — the 
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level fluctuations also change and are not of Wigner-Dyson type for 2k < m. 
We cannot pin down precisely the /c-value where such change occurs. The 
case m — 2, k — 1 suggests that this, too, happens near 2k = m. We have 
proved, however, that the level fluctuations become Poissonian in the limit 
k <C m <C I. Our analysis of the case k = 1 suggests that in this limit the 
eigenfunctions display localization in Fock space. With the analytical meth- 
ods at our disposal, we cannot reach that regime of the embedded ensembles 
where neither Poissonian nor Wigner-Dyson statistics applies. This is no 
surprise: We do not know of any case where such an aim would have been 
achieved analytically. 



It is interesting to compare our results with those of Ref. p9fl . There 
the two-point correlation function was calculated to lowest order in the loop 
expansion. (This expansion differs from the one in the present paper). It was 
found that the two-point function did have the dependence on r expected 
from canonical random matrix theory. Here, r is the difference in energy of 
the two Green functions in units of the mean level spacing. However, the 
coefficient multiplying r~ 2 turned out to be too large. It was speculated 
that higher-order terms in the loop expansion would reduce this coefficient 
to the value unity expected from canonical RMT. Now we know that this 
speculation applies only for 2k > m or so. In the limit k <C Tn -C I, the 
higher-order terms in the loop expansion of Ref. [2~9] must actually cancel 
the contribution of lowest order. 

Within the framework of the supersymmetry approach and in diagram- 
matic language, the transition at 2k = m is caused by the increased weight 
(with decreasing k) of intersecting Wick contraction lines. Universal random- 
matrix results are obtained whenever such contributions are negligible, and 
this apparently is the case for 2k > m. In other words, the intersecting 
contraction lines are responsible for the deviations from random-matrix uni- 
versality. We ascribe the special role played by the value 2k = m of the rank 
k to duality. The duality relation plays an important role in obtaining our 
results. It connects k with (m — k) and, thus, assigns a special role to the 
value 2k = m. 

From a physical point of view, our results can be understood in terms of 
the ratio Kp / P of the number of uncorrelated random variables over the num- 
ber of independent links, i.e., the number of non-zero entries in the matrix 
representation of the A;-body interaction in Hilbert space. We have shown 
that if all links were to carry uncorrelated random variables, the ensembles 
would have Wigner-Dyson type spectral fluctuations. Conversely, if all links 
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were to carry the same random variable, the ensembles would be completely 
integrable and display Poissonian statistics. These statements hold for all 
values of k. The actual situation is located between these two limits. The 
case k — 1 is closest to the integrable case, and the case k = m corresponds 
to canonical random-matrix theory. These facts show that deviations from 
random-matrix universality are not caused by the number of zeros in the 
matrix representation of the interaction but are strictly due to correlations 
between matrix elements. We may think of the random fc-body interaction 
as a sum of terms, each term carrying a single random variable multiplied by 
a fixed matrix. All random variables in the sum are uncorrelated. Reducing 
the number of independent variables is tantamount to regrouping this sum 
into a sum with a smaller number of terms. The question whether complete 
mixing between all states in Hilbert space is attained and random-matrix 
type results are obtained, is then a question of the number of terms in the 
sum. For instance, a stochastic coupling of only two fixed matrices will, in 
general, not lead to random-matrix results. 

In the light of these remarks, it seems appropriate to use the term local- 
ization for the results on eigenfunctions derived for k = 1 . Indeed, Anderson 
localization may be viewed as a problem of insufficient mixing in a stochastic 
Hamiltonian. This insufficient mixing is due to the restriction of the inter- 
action to nearest-neighbor couplings. If each lattice site were coupled to all 
other sites, localization would not occur. In the present context, it is not the 
limited number of links which causes localization but the further constraint 
on stochasticity caused by correlations between the matrix elements on these 
links. 

These results shed new light on the agreement of numerical results for 
EGOE(2) and EGUE(2), and of experimental data obtained in nuclear and 
atomic physics, with Wigner-Dyson spectral statistics. (We disregard here 
the complications which arise because the single-particle levels in major shells 
are non-degenerate in atoms and nuclei). Indeed, it is true both for the 
numerical work and for the experimental data that the underlying spaces are 
very far from the limit I — > oo considered in this paper. Among all these 
cases, the case of nuclei with excitation energies near neutron threshold is 
probably the one with the largest number of interacting configurations: The 
number of valence nucleons is around six or eight, and the total number of 
configurations is around 10 6 , implying values of I in the tens or twenties. A 
comparison of Figure |3| with Figure [1] shows that in this situation, the ratio 
Ki/P « 10~ 4 or so is relatively large. Thus, the fluctuations are near the 
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Figure 3: The ratio Ki/P on a logarithmic scale versus k for I = 25 and for 
various values of m. 

GOE or GUE limit although for I — > oo, none of these systems would display 
universal random-matrix behavior. Putting things very pointedly, we might 
say that the agreement between Wigner-Dyson level statistics and the results 
found in these systems, is a finite-size effect. 
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A Calculation of Loop Corrections 

A.l Loop Corrections for the EGUE(/c) 

The loop corrections are obtained in terms of an expansion around the 
saddle-point solution. The limit I — > oo is taken. For each value of (s, a), 
let denote the deviation of from the saddle-point value. First, 

we expand the logarithmic term in the effective Lagrangean in powers of 
5T, = J2 sa \^(k)C^5cr^ saS) and in powers of the source variable J. The 
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resulting terms are denoted by (p, q) where p indicates the power of 5S and 
q the power of J. Second, the resulting exponential is expanded again in 
powers of <5£ and J. We keep terms up to fourth order in <5£ and up to 
first (second) order in J for the one- (two)-point function, respectively. Fi- 
nally, we perform the Gaussian integration over 5a^ and identify among 
the remaining terms those which are of leading order in the limit I — > oo. 

One-point function: After expanding the logarithmic term, we retain 
the terms (p, q) with p = 1,...,4 and q = 0,1. The term (0,0) vanishes 
identically. Because of the saddle-point condition Eq. (f43|), the term (1,0) 
cancels the corresponding term from the expansion of (N/2) trg[cr^ a - ) ] 2 . 
After taking the trace over fi, we combine the term (2, 0) with the quadratic 
term from the expansion of (N/2) ^ sa trg[cr^ a ^] 2 and obtain a new Gaussian 
exponent of the form -(N/2) £Jl - (r' ') 2 (A^(£;)/A( )(A;))] trg [5a^] 2 . 
We are left with the terms (3,0), (4,0) and with the terms (p, 1) for all 
p — 0, ... ,4 as listed in Table |l| Expanding the exponential in powers of 
these terms, we keep terms up to fourth order in but omit terms odd in 

$(j( sa ) as these vanish upon the Gaussian integration. Moreover, we keep only 
terms linear in J and obtain [(0, 1) + (2, 1) + (4, 1) + (3, 0)(1, 1) + (4, 0)(0, 1)]. 
The term (4, 0)(0, 1) vanishes due to the Efetov-Wegner theorem. Evaluation 
of the remaining terms with the help of the orthogonality relations (p2|), 
the duality relation (|2T)|), and Wick's theorem for the Gaussian integrals is 
straightforward and yields non-zero contributions only from (0, 1) and from 
(4, 1). This gives the result in Eq. (HSI). 

Two-point function: After expanding the logarithm, we retain all terms 
(p,q) with p = 1,...,4 and q = 0,1,2. Again, (0,0) vanishes identically, 
(1, 0) cancels and (2, 0) forms part of the new Gaussian exponent. For con- 
venience, we keep the terms (0, 1) and (0, 2) in the exponent as well. We 
are left with the terms (3,0), (4,0), (p, 1) and (p, 2) for p = 1,...,4 as 
shown in Table 0. The exponential is expanded in a power series where we 
keep again only terms even and up to fourth order in ££, and terms up 
to second order in the source variable J. Performing the Gaussian integra- 
tion and dropping terms which vanish as 1/N for iV — > oo, we are finally 
left only with the contributions stemming from (4,1) and (2,1)(2,1). We 
rewrite these terms as an exponential. From (4, 1), we obtain a correction to 
the canonical random-matrix contribution as given in expression (0). The 
only non-trivial correction stems from (2, 1)(2, 1) and is proportional to the 
expression ([52]). 
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p q 



3 -§ (^S) 3 tr,trg(5£) 3 

4 -I (g}) 4 tr M trg(5S) 4 

1 ~xSt tr^trgj 

1 1 -(w>) 2 tr M trg(5SJ) 

2 1 -(SJ) 3 tr,trg((5E) 2 J) 

3 1 -(S) 4 tr,trg((5S) 3 J) 

4 1 -(S) 5 tr M trg((5S) 4 J) 

Table 1: Relevant terms in the expansion of the logarithmic term in the 
effective Lagrangean of the EGUE(A;) one-point function. Parameters p and 
q indicate powers of 5T, and J, respectively. 

A. 2 Loop Corrections for the Limiting Ensemble 

EGUE max (/c) 

For the limiting ensemble EGUE max (A;), the loop corrections are calculated 
in very much the same way as for the EGUE(A;). It is important, however, to 
keep in mind the operator character of a in the Hilbert space of the (m — k)- 
particle system. The logarithmic term in the effective Lagrangean contains 
the matrix elements (fj,\A^SaA\fj) of 5a. The calculation is carried through 
independently for each value of //. After expansion of the logarithm, it is 
convenient to include in the expansion of the exponent also the term (2,0). 
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P Q 



3 
4 
1 

2 
3 
4 



2 
3 
4 





1 

1 



!(A(°))~ 3 tr M trg(L 5£) 3 
i(A(°))" 4 tr M trg(L 5£) 4 
(A^)- 2 tr M trg(L 5E L J) 

i(\W)- 3 tr M trg f(L 5£) 2 L J 



1 -(A(°))- 4 tr M trg ((L 5E) 3 L j) 
1 -i(A(°))- 5 tr^trg ((L L j) 



2 i(A (0) )~ 3 tr M trg (L (L J) 



(A(°))- 4 tr M trg (L 5£) 2 (L J) 2 + ±(L L J) ; 
^(A(°))- 5 tr^trg \(L 5S) 3 (L J) 2 + (L L J) 2 



2 (A(°))- 6 tr M trg (L 5S) 4 (L J) 2 + (L 5S) 2 (L L J) 2 



+ 



±((L 5S) 2 L j)' 



Table 2: Relevant terms in the expansion of the logarithmic term in the 
effective Lagrangean of the EGUE(/c) two-point function. Parameters p and 
q indicate powers of 5Y, and J, respectively. The expressions are evaluated 
at the center of the spectrum, E = 0. We put J = (e/2)TLT _1 + TJT~\ 
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We again restrict ourselves to terms up to fourth order in 5a and up to first 
or second order in J. 

One-point function: After expansion of the exponential, we are left with 
[(0, 1) + (2, 1) + (4, 1) + (2, 0)(0, 1) + (2, 0)(2, 1) + (3, 0)(1, 1) + (4, 0)(0, 1)]. 
The terms (2,0) (0,1) and (4,0) (0,1) yield zero due to the Efetov-Wegner 
theorem. After the Gaussian integration we are left with contributions from 
the canonical term (0, 1) and from (4, 1). Both together yield 

■ <62) 

The correction to the leading term is ((™) / 'N {k)) A : . Even in the worst case, 
k = 1, this contribution vanishes like /~ 4 for I — > oo. 

Two-point function: After expanding the exponential and performing 
the Gaussian integration, we are left with the terms [(0, 2) + (2, 2) + (4, 2) + 
(2, 0)(2, 2) + (0, 1)(0, 1) + (0, 1)(4, 1) + (1, 1)(1, 1) + (1, 1)(3, 1) + (2, 1)(2, 1) + 
(2, 0)(1, 1)(1, 1)]. Except for the canonical contribution (0, 1)(0, 1), all these 
terms vanish at least like ((™)/A(A;)) 4 , as does the correction term for the 
one-point function. 



Nt (0) 

■L " ' max 

N(k) 



trgj 



B Binary Correlation Approximation for the 
n— Point Funtion 

In this Appendix, we show that for k <^ m ^ I, the eigenvalues of EGUE(A;) 
have a Poissonian distribution. The proof applies analogously to EGOE(/c). 
Our proof extends the calculation of the average one-point function of Section 
5.2 of Ref. |29| to the case of the n-point function with n > 1. In order to 
avoid the occurrence of terms which diverge in the limit / — > oo in the 
sums introduced below, we choose in this Appendix the second moment of 
the matrix elements of V& differently from the main body of the paper (cf. 
Eq. (||)) by putting Vq = (A^/c)) -1 . This choice has the merit of yielding 
an average Gaussian spectrum with unit width for all choices of k, m, and I. 
We briefly recall the method of Ref. |[29| . The average one-point function 
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is expanded in inverse powers of z and the summation and the average over 
the ensemble are interchanged, 

i 00 1 

^) = iv^^ tr -(^)- (64) 

p=0 

Because of the Gaussian distribution of 14, only even powers of V& contribute. 
The average is performed via Wick contraction. There are (2p — 1)!! ways of 
pairwise contracting the V^'s. For k <^ m <^ I, each contracted pair of \4's 
contributes a factor unity. The result is 

1 00 1 

^) = ^E^T^- 1 ) !! - ( 65 ) 

p=0 

Borel summation of this series gives a closed expression for g(z) and the 
Gaussian form for the average level density, see Ref . . 

In considering next the average two-point function g(zi)g(z 2 ), we use the 
expansion ( |B"4T) for g(zi) and for g(z 2 ). With z\ = E + +e/2 and z 2 = E~—e/2, 
the letters E and e denote the mean value and the difference of the energy 
arguments of the two Green functions, respectively. An upper plus (minus) 
sign denotes an infinitesimal positive (negative) imaginary energy increment. 
We are interested in values of e which are of the order of the mean level 
spacing and, thus, small compared to \E\ with \E\ ^> 1. (We recall that the 
spectrum has width unity). Then, we have approximately \z\\ ~ \z 2 \. In the 
expansion, we collect all terms which carry the same even power 2n of V}. . 
This yields 

gizMz,) =EE 7 m 7W ^IW'] tr u [(V k )^] . (66) 

„=0 t=0 Z l Z 2 ly 

Wick contraction may lead to s = 0, 1, . . . , n cross-contracted pairs of V^'s, 
i.e., pairs where the two V^'s reside in different traces. The remaining number 
of \4's in each trace must be even. Hence, 

00 n n—s ^ 

s+2p+l 2n+l-s-2p 
n=0 s =0 p=0 Z l Z 2 

x ^(tr M [(^) s+2p ] tr4(\4) 2 «— 2p]^ . (67) 
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The index s indicates that s pairs are cross-contracted. The number of 
different ways of cross-contracting the V^'s is ( s ~^ 2p ) ( 2n ~f'~ s ), while there are 
(2p— 1)!! (2n — 2p— 1)!! ways of pairwise contracting the remaining V^s. As 
in the case of the one-point function, each of the latter contractions yields 
the factor unity, irrespective of the sequence in which the V&'s appear under 
the trace, provided that C m < L Likewise, for k <C m <C I each way 
of cross-contracting the remaining s pairs of 14 's yields the same answer, 
irrespective of where they appear in each of the traces. The result is 

— -— - ^ " — C+2P) (2n-2p- a) (2p _ 1)n(2w 2p _ 1)H 

g\Z 1 )g{Z 2 ) - 2-^2.^1^ s+2p+l 2n+l-s-2p 

n=0 s=0 p=0 Z l Z 1 

x ^ (t#yi«l) s • (68) 



We now show that for each term with fixed n in Eq. (pSj), in the limit 
/ — > oo the terms with s ^ become negligibly small in comparison with 
the terms that have s = 0. The latter correspond to the unlinked contribu- 
tions. Upon resummation, these terms yield the product g(zi) ■ g(z 2 ). This 
shows that the two-point function factorizes, and that the linked correlation 
function R 2 (zx,z 2 ) defined in Eq. ([H]) vanishes asymptotically in the limit 
/ — > 00. 

To this end, we observe that the binomial factors and powers of z\ and z 2 
in Eq. fl68|) do not depend on I or on m. Moreover, we recall that \z\\ ~ | ^2 1 - 
Hence, the only /-dependent part in the sum in Eq. (pE|) is 



T ' = Jpfa[<yk)']te»[(V k )']) M . (69) 

Obviously, we have T = 1. Thus, T and with it all unlinked contributions 
are constant and independent of /. To prove our claim, it suffices to show 
that T s vanishes at least as l~ l for all s > 1. 

We write 14 = Yliap v apA J[ a Ap ) with v a p and A\ a short-hand notation for 
the random variables and for a] 1 • • • dj k with j \ < j 2 < . . . < jk, respectively, 
and likewise for Ap. Then, 



Ts ~ N 2 (A^(k)) s tx n( A U A Pi ■ ■ ■ A^ApJtr^A^A^ . . . A\A C 



ai,f3i,...,as,/3s 



(70) 
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Without loss of generality, the indices a.\\ a 2 ; ■ ■ ■ ; a s can be chosen as 1, . . . , k; 
k + 1, . . . , 2k; . . . ; (s — l)k + 1, . . . , sk. This yields an overall factor . We 



find 

'i 



( ) 

Ts = N 2 (A(°)(k)) s ^ tT v( A \,...,k A Pi ■ ■ ■ A \ s -i)k+i,...,sk A Ps) 
x tx v {A\ i A 1 ^ k ...A\A {s _ 

l)k+l,...,sk) ■ (71) 

The summation over the indices j3i, . . . , /3 S is now restricted to a sum over 
all permutations of the integers 1,2, ... ,sk. We note that the number (sk)\ 
of such permutations does not depend on / or m. We consider one such 
permutation and show that the resulting dependence on I is the same for all 
permutations. 

In both traces, any permutation generates terms of the type a^a, and 
terms of the type ajOj. Let the number of terms of the first type in the 
first trace be i, those of the second type in the first trace be sk — i. It is 
easy to see that the corresponding numbers in the second trace are sk — i 
and i, respectively. This implies that in the first trace, only those states 
contribute for which in the first sk single-particle states, i states are occupied. 
The remaining m — i particles must be distributed over the remaining / — sk 
single-particle states. This yields a total of Hilbert-space states which 

contribute to the first trace. The corresponding counting for the second trace 
yields ( m Z _J s s fc fc +i ). Altogether, we find that the contribution of this particular 
permutation to the quantity T s is given by 

HY (l-sk\( l-sk j ^ 



7V 2 (A(°)(fc)) s \m - ij \m - sk + i 

Using Stirling's formula, we find that for / — > oo, this expression behaves as 
l~ sk , independent of the value of i. This then shows that all permutations 
yield the same behavior. Their number is independent of / and m. (Actually, 
the single-particle operators appear in different sequences, giving perhaps 
rise to cancellations. This possibility is not considered here). We conclude 
that all connected terms in the two-point function vanish at least as fast 
as l~ k , so that for / — » oo, the two-point function is given by the product 
of the averaged one-point functions and the two-point correlation function 
^2(^1,^2) in Eq. (|6l|) vanishes. 

The same line of reasoning can be applied to the n-point functions with 
n > 2. We find that every totally connected piece in the n-point function 
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containing s cross-contracted pairs of V^'s tends to zero as l~ sk for I — > oo. 
This shows that all higher n-point functions are asymptotically likewise given 
by the product of n average Green functions. Thus, all correlation functions 
vanish asymptotically, and the spectral fluctuations for EGOE(fc) become 
Poissonian in the limit k <m <C Z. 
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